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In  this  paper,  a capacitated  vehicle  routing  problem  is  discussed  which  occurs  in  the  context  of  glass  waste 
collection.  Supplies  of  several  different  product  types  (glass  of  different  colors)  are  available  at  customer 
locations.  The  supplies  have  to  be  picked  up  at  their  locations  and  moved  to  a central  depot  at  minimum 
cost.  Different  product  types  may  be  transported  on  the  same  vehicle,  however,  while  being  transported 
they  must  not  be  mixed.  Technically  this  is  enabled  by  a specific  device,  which  allows  for  separating  the 
capacity  of  each  vehicle  individually  into  a limited  number  of  compartments  where  each  compartment  can 
accommodate  one  or  several  supplies  of  the  same  product  type.  For  this  problem,  a model  formulation  and 
a variable  neighborhood  search  algorithm  for  its  solution  are  presented.  The  performance  of  the  proposed 
heuristic  is  evaluated  by  means  of  extensive  numerical  experiments.  Furthermore,  the  economic  benefits  of 
introducing  compartments  on  the  vehicles  are  investigated. 

© 2015  Elsevier  B.V.  and  Association  of  European  Operational  Research  Societies  (EURO)  within  the 
International  Federation  of  Operational  Research  Societies  (IFORS).  All  rights  reserved. 


1.  Introduction 

The  vehicle  routing  problem,  which  will  be  discussed  in  this  paper, 
is  a variant  of  the  classic  capacitated  vehicle  routing  problem  (CVRP; 
for  surveys  see  Golden,  Raghavan,  & Wasil,  2008;  Laporte,  2009;  or 
Toth  & Vigo,  2014)  and  occurs  in  the  context  of  glass  waste  collection 
in  Germany.  Glass  waste  has  to  be  recycled  by  law  and  is  used  as  a raw 
material  for  the  production  of  new  glass  products.  It  has  to  be  taken  to 
recycling  stations  by  the  consumers  where  it  is  disposed  into  differ- 
ent containers  according  to  the  color  of  the  waste  (usually  colorless, 
green  and  brown  glass).  Colors  are  kept  separated  because  the  pro- 
duction of  new  glass  products  is  less  cost-intensive  if  the  glass  waste 
is  not  too  inhomogeneous  with  respect  to  its  color.  Trucks,  which  are 
located  at  a depot  of  a recycling  company,  pick  up  the  glass  waste 
from  the  recycling  stations.  Since  they  possess  a relatively  large  load- 
ing capacity,  they  can  call  at  several  recycling  stations  before  they 
have  to  return  to  the  depot.  Recent  truck  models  are  equipped  with  a 
special  device  which  allows  for  introducing  bulkheads  in  predefined 
positions  of  the  loading  space  such  that  it  can  be  split  into  different 
compartments  and,  thus,  enabling  transportation  of  glass  waste  with 
different  colors  on  the  same  truck  without  mixing  the  colors  on  the 
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tour.  This  gives  rise  to  the  question  how  the  tours  of  the  trucks  should 
be  designed  given  the  availability  of  a device  of  this  kind. 

The  problem  under  discussion  can  be  classified  as  a multi- 
compartment vehicle  routing  problem  (MCVRP).  However,  it  is  differ- 
ent from  the  ones  previously  discussed  in  the  literature  with  respect 
to  the  following  properties: 

• The  size  of  each  compartment  is  not  fixed  in  advance  but  can  be 
determined  individually  for  each  vehicle/each  tour. 

• The  size  of  the  compartments  can  only  be  varied  discretely,  i.e. 
the  walls  separating  the  compartments  from  each  other  can  only 
be  introduced  in  specific,  predefined  positions. 

• The  number  of  compartments,  into  which  the  capacity  of  a vehicle 
is  divided,  can  be  identical  to  the  number  of  product  types  (glass 
waste  types)  but  can  also  be  smaller. 

Consequently,  not  only  the  vehicle  tours  have  to  be  determined,  but 
it  has  also  to  be  decided  for  each  vehicle/tour  (i)  into  how  many  com- 
partments the  vehicle  capacity  should  be  divided,  (ii)  what  the  size 
of  each  compartment  should  be,  and  (iii)  which  product  type  should 
be  assigned  to  each  compartment. 

The  problem  is  NP-hard,  since  it  is  a generalization  of  the  CVRP 
(see,  for  example,  Toth  & Vigo,  2014).  By  application  of  a mathe- 
matical model-based  exact  solution  approach,  we  were  only  able 
to  solve  problem  instances  with  a limited  size  to  optimality.  There- 
fore, a heuristic,  namely  a variable  neighborhood  search  (VNS),  has 
been  developed  and  will  be  presented.  According  to  the  best  of  our 
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knowledge,  this  is  the  first  method  which  has  been  proposed  for  this 
problem  so  far.  We  will  further  analyze  what  the  economic  benefits 
are  which  stem  from  the  introduction  of  flexibly  sizable  compart- 
ments. 

The  remainder  of  this  paper  is  organized  as  follows.  Section  2 
presents  a formal  definition  and  a mathematical  formulation  of  the 
problem.  The  relevant  literature  related  to  the  MCVRP  is  discussed  in 
Section  3.  In  Section  4,  the  proposed  variable  neighborhood  search 
algorithm  is  introduced.  Extensive  numerical  experiments  have  been 
performed  in  order  to  evaluate  the  mathematical  model  and  the  VNS. 
The  design  of  these  experiments  and  the  corresponding  results  are 
presented  in  Section  5.  Finally,  the  main  findings  are  summarized  and 
an  outlook  on  future  research  is  given  in  Section  6. 

2.  Problem  description  and  formulation 

The  multi-compartment  vehicle  routing  problem  with  flexible 
compartment  sizes  (MCVRP-FCS)  can  be  formulated  as  follows:  Let 
an  undirected,  weighted  graph  G = (V,  E ) be  given  which  consists  of 
a vertex  set  V = {0, 1, . . . , n},  representing  the  location  of  the  depot 
({0})  and  the  locations  of  n customers  ({1, . . . , n}),  and  an  edge  set 
E = { (i,  j):i,  j e V,  i < j},  representing  the  edges  which  can  be  trav- 
eled between  the  different  locations.  To  each  of  these  edges,  a non- 
negative cost  Cfj,  ( i , j)  g E,  is  assigned.  It  is  assumed  that  all  of  these 
costs  satisfy  the  triangle  inequality. 

Further,  let  a set  P of  product  types  be  given.  At  each  vertex  (ex- 
cept for  the  depot)  exists  a non-negative  supply  sip(i  e V\{0},  p e P) 
of  each  of  the  product  types.  The  supplies  have  to  be  collected  at  their 
locations  and  transported  to  the  depot  without  the  product  types  be- 
ing mixed.  A location  may  be  visited  several  times  in  order  to  pick 
up  different  product  types.  However,  if  being  picked  up,  each  supply 
has  to  be  loaded  in  total.  In  other  words,  a split  collection  of  a single 
supply  is  not  permitted. 

For  the  purpose  of  transportation,  a set  I<  of  homogeneous  ve- 
hicles is  available,  each  equipped  with  a total  capacity  Q.  Individu- 
ally for  each  vehicle  k e I<,  the  total  capacity  Q.  can  be  divided  into 
a limited  number  m of  compartments,  m < \P\,  which  allows  for 
loading  products  of  different  types  on  a single  vehicle  while  keep- 
ing them  separated  during  transportation.  The  size  of  the  compart- 
ments can  be  varied  discretely  in  equal  step  sizes,  i.e.  each  com- 
partment size,  but  also  the  total  vehicle  capacity  Q,  is  an  integer 
multiple  of  a basic  compartment  unit  size  qumt.  Let  the  set  of  these 
multiples  be  denoted  by  M = {0,  1,2,  . . . , mmax}  where  mmax  = 
QJqumt.  Then  qm  = ^ax  • m,  m e M,  denotes  a compartment  size 
relative  to  the  total  capacity  Q.  which  consists  of  m (m  e M ) multi- 
ples of  the  basic  compartment  unit  size  qunit.  To  illustrate  these  as- 
pects, we  introduce  a small  example  in  which  the  vehicle  capacity  Q. 
amounts  to  200  units  and  the  basic  compartment  unit  size  qunit  to  10 
units.  Hence,  only  compartment  sizes  of  10,  20,  30,...,  200  units  or 
5 percent,  10  percent,  15  percent,...,  100  percent  of  the  vehicle  ca- 
pacity can  be  selected.  Accordingly,  mmax  is  equal  to  20  and  a com- 
partment with  m = 7 corresponds  to  a relative  compartment  size  of 
Qm  = m max  • m = jq  • 7 = 0.35,  i.e.  35  percent  of  the  vehicle  capac- 
ity. It  is  important  to  note  that  the  set  of  potential  compartment  con- 
figurations is  identical  for  all  vehicles.  However,  the  actual  configura- 
tion in  a particular  solution  might  be  different  for  each  vehicle. 

What  has  to  be  determined  is  a set  of  vehicle  tours,  an  assignment 
of  product  types  to  the  vehicles  and  the  sizes  of  the  corresponding 
compartments  such  that  all  supplies  are  collected,  that  the  capacity 
of  none  of  the  used  vehicles  is  exceeded,  and  that  the  total  cost  of  all 
edges  to  be  traveled  is  minimized. 

This  problem  involves  the  following  partial  decisions  to  be  made 
simultaneously: 

• assignment  of  product  types  to  each  of  the  vehicles 
(this  decision  determines  which  product  types  can  be  collected  by 
each  vehicle); 


• determination  of  the  size  of  each  compartment 

(this  decision  fixes  for  each  vehicle  how  its  total  capacity  is  split 
into  compartments); 

• assignment  of  supplies  to  each  of  the  vehicles 

(this  decision  implicitly  includes  an  assignment  of  locations  to 
vehicles); 

• sequencing  of  the  locations  for  each  of  the  vehicles 

(this  decision  determines  for  each  vehicle  in  which  sequence  the 
assigned  locations  are  to  be  visited). 

We  note  that  every  vehicle  routing  problem  involves  decisions 
of  the  last  two  types,  while  the  first  and  the  second  one  define  the 
uniqueness  of  the  MCVRP-FCS. 

In  order  to  formulate  a mathematical  model  for  the  MCVRP-FCS, 
we  introduce  the  following  four  types  of  variables: 


if  supply  of  product  type  p at  location  i is  collected 
by  vehicle  k, 
otherwise, 

ie  V\{0},peP,keI<; 

if  i = 0 and  edge  (i,  j)  is  used  twice  by  vehicle  k, 
if  edge  (i,  j)  is  used  once  by  vehicle  k, 
otherwise, 

i,  j eV  :i<  j,  k e I<; 

if  size  qm  is  selected  for  product  type  p in 
vehicle  k, 
otherwise, 

peP,keK,me  M; 

if  location  i is  visited  by  vehicle  k, 
otherwise, 

i eV,ke  I(. 


The  objective  function  and  the  constraints  of  the  model  can  then 
be  formulated  as  follows: 

min  E E CijXiJk  (1) 

(i,j)(zA  keK 


E U‘Pk  = 1 Vi’  6 'AW’  P 6 P : SW  > 0 

keK 

(2) 

Uipk  < zik  Vi  e V\{0},  P 6 p,  k e K 

(3) 

Zik  < Z0k  Vi  € V\{0},  k € K 

(4) 

J2JJx0)k<2\I(\ 

j^v-  keK 
i>  o 

(5) 

E + E xi ik  = 2z"<  WcV’keK 

jeV:  jeV: 

i<  j j<i 

(6) 

E E ypkm  < ™ W 6 K 

peP  meM 

(7) 

E E <3  ™ypkm  <1  Vk  € K 

peP  meM 

(8) 

E sipuiPk<Q-J2Qmypkm  VpeP,keI<  (9) 

ieV\{0}  meM 


EE*y*  < |S|  - 1 Vk  e K,S  c V\{0}  : |S|  > 2 (10) 

ieS  J'eS: 
i<j 
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a)  A MCVRP-FCS  instance  with  three 
locations  and  four  product  types 


b)  A solution  of  the  MCVRP-FCS  instance 


Fig.  1.  Illustrative  example  of  a Multi-Compartment  Vehicle  Routing  Problem  with  Flexible  Compartment  Sizes. 


uipke{  0,1}  WieV\{0],peP,keI<  (11) 

xojk  e {0, 1,  2}  Vj  e V\{0},  k e I<  (12) 

xijk  e {0, 1}  Vi  g V\{0},  j e V\{0}  : i < j,  k e I(  (13) 

yPkm£  {0,1}  VpeP,keK,meM  (14) 

zike{  0,1}  VieV,keK  (15) 


The  objective  function  (1 ) defines  the  total  cost  of  all  tours,  which 
has  to  be  minimized.  Constraints  (2)  guarantee  that  each  positive 
supply  is  assigned  to  exactly  one  vehicle,  while  constraints  (3)  ensure 
that  a supply  may  only  be  satisfied  by  a vehicle  if  the  corresponding 
customer  location  is  visited  by  this  vehicle.  Constraints  (4)  make  sure 
that  the  depot  is  included  in  each  tour.  According  to  constraints  (5), 
the  number  of  tours  is  limited  by  the  number  of  available  vehicles. 
Constraints  (6)  represent  the  vehicle  flow  constraints.  Constraints  (7) 
ensure  that  the  number  of  product  types  assigned  to  a vehicle  may 
not  exceed  the  maximum  number  of  compartments.  Constraints  (8) 
ensure  that  for  each  vehicle  the  sum  of  the  relative  compartment 
sizes  qm  must  be  smaller  than  or  equal  to  1 (100  percent),  i.e.  the  sum 
of  the  actual  compartment  sizes  must  be  smaller  than  or  equal  to  the 
vehicle  capacity.  In  addition  to  that,  constraints  (9)  guarantee  for  each 
compartment  that  its  capacity  is  not  exceeded  by  the  assigned  sup- 
plies. Constraints  (10)  represent  the  subtour-elimination  constraints. 
Finally,  constraints  (11 )— (15)  characterize  the  variable  domains. 

An  illustrative  example  for  the  MCVRP-FCS  with  three  locations 
(vlt  v2,  v3)  and  four  product  types  (R:  red,  Y:  yellow,  G:  green,  B:  blue) 
is  provided  in  Fig.  la.  The  locations  are  depicted  as  circles  while  the 
product  types  are  depicted  as  squares.  The  number  within  a square 
indicates  the  amount  of  supply  of  a certain  product  type  at  a certain 
location.  In  case  no  supply  exists,  the  respective  square  is  crossed.  In 
the  shown  example,  each  location  has  a supply  of  one  unit  for  each 
product  type  with  the  exception  of  location  3 not  having  any  supply 
of  product  type  B at  all.  Furthermore,  each  vehicle  has  a capacity  of 
six  units  and  the  maximum  number  of  compartments  equals  three. 

A solution  to  this  problem  instance  with  two  vehicles  is  shown 
in  Fig.  lb.  For  each  vehicle,  a distinguished  panel  is  depicted,  which 
consists  of  information  about  the  product-vehicle-assignment,  the 
supply-vehicle-assignment,  and  the  location  sequencing.  If  a certain 


product  type  is  assigned  to  a vehicle,  the  associated  square  in  the 
top  left  part  of  the  panel  is  colored  and  named  accordingly.  Further, 
if  a certain  supply  is  assigned  to  a vehicle,  the  corresponding  square 
below  the  related  location  is  colored  and  includes  the  amount  of 
the  supply.  Finally,  the  routes  are  indicated  by  edges  between  the 
vertices. 

3.  Literature  review 

Some  variants  of  the  MCVRP  have  been  discussed  in  the  literature 
already.  El  Fallahi,  Prins  and  Wolfer  Calvo  (2008)  describe  a MCVRP 
arising  in  the  context  of  distribution  of  animal  food  to  farms  where 
the  sizes  of  the  compartments  and  the  corresponding  assignment  of 
product  types  to  the  compartments  are  fixed  in  advance.  They  pro- 
pose a memetic  algorithm  and  a tabu  search  algorithm  for  solving 
this  problem.  Muyldermans  and  Pang  (2010)  consider  a similar  prob- 
lem and  develop  a guided  local  search  algorithm.  Furthermore,  they 
identify  problem  characteristics  under  which  transportation  of  sev- 
eral product  types  on  vehicles  with  multiple  (a  priori  fixed)  compart- 
ments outperforms  separate  transportation  of  the  different  product 
types  on  vehicles  with  a single  compartment.  Chajakis  and  Guignard 
(2003)  present  a similar  problem  in  the  context  of  deliveries  to  con- 
venience stores.  For  their  problem,  they  develop  heuristics  based  on 
Lagrangian  Relaxations.  In  contrast  to  these  contributions,  in  this  pa- 
per a MCVRP  with  flexible  compartment  sizes  and  no  a priori  assign- 
ment of  product  types  to  compartments  is  considered. 

Avella,  Boccia  and  Sforza  (2004)  and  Brown  and  Graves  (1981 ) de- 
scribe applications  of  the  MCVRP  in  the  context  of  petrol  replenish- 
ment. In  these  problems,  the  compartment  sizes  are  fixed  in  advance, 
while  the  assignment  of  product  types  to  compartments  is  not.  In 
contrast  to  the  MCVRP-FCS,  only  one  supply  can  be  assigned  to  each 
compartment,  because  petrol  tanks  must  be  emptied  completely  dur- 
ing the  visit  of  one  customer.  Moreover,  Cornillier  et  al.  (2008)  extend 
the  problem  to  multiple  periods  with  inventory  components. 

Derigs  et  al.  (2011)  introduce  and  examine  two  variants  of  the 
MCVRP:  The  first  variant  is,  in  general,  similar  to  the  problem  de- 
scribed by  El  Fallahi  et  al.  (2008),  Muyldermans  and  Pang  (2010),  and 
Chajakis  and  Guignard  (2003),  whereas  the  second  variant  is  charac- 
terized by  continuous  flexible  compartment  sizes,  while  the  assign- 
ment of  product  types  to  compartments  is  not  fixed.  However,  in  con- 
trast to  the  MCVRP-FCS,  each  customer  may  have  multiple  demands 
for  the  same  product  type.  The  case  in  which  the  number  of  com- 
partments may  be  smaller  than  the  number  of  product  types  is  not 
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explicitly  considered.  The  authors  introduce  a general  mathematical 
model,  i.e.  one  which  covers  both  problem  variants,  and  a large  neigh- 
borhood search  algorithm  which  solves  their  variants  of  the  MCVRP. 

A problem  similar  to  the  MCVRP-FCS  is  presented  by  Caramia  and 
Guerriero  (2010).  They  present  a case  study-oriented  MCVRP  with 
heterogeneous  vehicles  in  milk  collection  with  fixed  compartment 
sizes  but  no  a priori  assignment  of  product  types  to  vehicles.  Also, 
the  number  of  product  types  (four  in  the  case  study)  may  be  smaller 
than  the  number  of  compartments  per  vehicle  (three  or  more  com- 
partments). In  contrast  to  the  MCVRP-FCS,  they  also  consider  trailers 
which  can  be  added  to  a vehicle.  For  their  problem,  they  present  a 
model-oriented  heuristic  procedure  which  assigns  demands  to  vehi- 
cles in  a first  step  and  determines  specific  routes  for  these  vehicles  in 
a second  step.  Only  recently,  Archetti,  Campbell  and  Speranza  (2014) 
presented  theoretical  comparisons  between  multi-commodity  and 
single-commodity  vehicle  routing  problems,  where  they  provide  the- 
oretical and  empirical  insights  on  relationships  between  the  split  de- 
livery vehicle  routing  problem,  the  multi-compartment  vehicle  rout- 
ing problem,  and  a combined  split  delivery  and  multi-compartment 
vehicle  routing  problem.  Finally,  Mendoza  et  al.  (2010,  2011)  discuss 
the  multi-compartment  vehicle  routing  problem  with  stochastic  de- 
mands for  which  they  propose  several  construction  procedures  and  a 
memetic  algorithm. 

Apart  from  the  mentioned  applications  in  glass  waste  collection, 
fuel  delivery,  and  milk  collection,  multi-compartment  routing  prob- 
lems are  found  in  several  other  contexts.  El  Fallahi  et  al.  (2008)  de- 
scribe a problem  occurring  in  the  distribution  of  animal  food  to  farms, 
in  which  the  food  for  a certain  species  should  always  be  assigned 
to  the  same  compartments  because  of  sanitary  rules.  Chajakis  and 
Guignard  (2003)  deal  with  the  simultaneous  delivery  of  dry,  refriger- 
ated, and  frozen  products  to  convenience  stores.  Lahyani  et  al.  (2015) 
present  a case  study  from  Tunisia  in  which  olive  oil,  categorized  into 
different  types  of  quality  grades,  has  to  be  collected  from  produc- 
ers in  separate  compartments.  Finally,  multi-compartment  problems 
also  occur  in  the  context  maritime  transportation  when  bulk  prod- 
ucts have  to  be  transported  (Fagerholt  & Christiansen  2000). 

Our  paper  attends  to  additional  aspects  which  - according  to  the 
best  of  our  knowledge  - have  not  been  considered  in  detail  in  the 
literature  before.  These  aspects  include  the  case  of  discretely  flexible 
compartment  sizes  and  the  case  in  which  the  number  of  compart- 
ments per  vehicle  is  smaller  than  the  number  of  products  types.  Es- 
pecially the  second  aspect  defines  a new,  unique  problem  structure 
for  which  we  suggest  a heuristic  solution  approach. 

4.  Variable  neighborhood  search 

4.1  Overview 

In  order  to  determine  good  solutions  for  large  problem  instances 
of  the  MCVRP-FCS,  a variable  neighborhood  search  approach  with 
multiple  starts  (MS-VNS)  has  been  developed.  The  concept  of  VNS 
was  first  proposed  by  Mladenovic  and  Hansen  (1997)  and  Hansen 
and  Mladenovic  (2001).  VNS  is  a local  search-based  metaheuristic 
which  explores  the  solution  space  by  means  of  multiple  neighbor- 
hood structures. 

The  general  VNS  framework  can  be  described  as  follows:  Given 
a set  of  different  neighborhood  structures  which  are  sequenced  in  a 
specific  order,  VNS  starts  from  an  initial  solution  x (also  first  incum- 
bent solution)  and  randomly  selects  a neighbor  x'  of  x according  to  the 
first  neighborhood  structure  (shaking).  The  neighbor  is  subsequently 
improved  by  application  of  a local  search  algorithm,  providing  a solu- 
tion x"  (improvement  phase  #2 ).  If  this  solution  does  not  satisfy  a given 
acceptance  criterion,  another  neighbor  x!"  of  the  incumbent  x solu- 
tion is  selected  randomly,  this  time,  however,  from  the  neighborhood 
defined  by  the  next  neighborhood  structure  in  the  sequence,  and  the 
above-described  process  is  repeated.  Otherwise,  if  x!'  satisfies  the  ac- 


ceptance criterion,  it  is  updated  as  the  new  incumbent  solution,  i.e. 
x :=  x!\  and  the  index  of  the  neighborhood  structure  is  set  back  to  1, 
i.e.  for  the  next  pass  of  the  improvement  phase  #1,  a neighbor  x’  of  the 
new  incumbent  solution  of  x is  selected  from  the  neighborhood  de- 
fined by  the  first  neighborhood  structure.  Should  no  new  incumbent 
solution  be  determined  over  a run  through  all  neighborhood  struc- 
tures, the  procedure  restarts  from  the  first  neighborhood  structure  in 
the  sequence.  The  procedure  continues  until  a predefined  termina- 
tion criterion  is  attained.  In  the  following,  the  execution  of  a shaking 
step  and  the  subsequent  improvement  phase  will  be  referred  to  as 
one  iteration. 

Usually  the  sequence  of  the  neighborhood  structures  is  chosen  in  a 
way  that  the  size  of  the  neighborhoods  increases  as  the  VNS  proceeds. 
As  a consequence,  initially  the  procedure  explores  neighbors  from  the 
incumbent  solution,  which  can  be  obtained  by  small  changes,  while 
it  only  moves  on  to  explore  larger  changes  if  the  search  for  better 
solutions  was  unsuccessful  in  such  “close”  neighborhoods. 

As  we  have  mentioned  before,  an  important  aspect  of  the  MCVRP- 
FCS  are  the  four  different  types  of  partial  decisions  which  have  to  be 
made  (see  Section  2).  The  usage  of  a single  neighborhood  structure 
may  revise  several  of  these  decisions  simultaneously.  However,  the 
combination  of  decisions  is  often  identical.  By  the  usage  of  more  than 
one  neighborhood  structure,  VNS  is  able  to  revise  different  combina- 
tions of  problem  decisions,  which  permits  a more  diverse  exploration 
of  the  solution  space. 

One  of  the  disadvantages  of  this  algorithm  is  that  it  often  tends 
to  get  stuck  in  a local  optimum.  In  order  to  overcome  this  drawback, 
we  implemented  our  VNS  in  connection  with  a multi-start  approach. 
The  VNS  restarts  several  times  from  different  initial  solutions,  which 
are  generated  by  a randomized  construction  procedure.  The  num- 
ber of  starts  (in  the  following  called  loops)  is  also  limited  by  a ter- 
mination criterion  (ms  termination  criterion).  Furthermore,  each  time 
a new  best  solution  has  been  found,  it  is  attempted  to  improve  this 
solution  by  a second  improvement  phase  (improvement  phase  #2). 

A pseudo-code  of  our  VNS  procedure  is  presented  in  Fig.  2.  Zbest 
denotes  the  objective  function  value  of  the  best  solution  found  over 
all  loops  and  Z(x)  denotes  the  objective  function  value  of  an  arbitrary 
solution  x.  Details  of  the  randomized  construction  procedure,  the  so- 
lution space,  the  neighborhood  structures  and  improvement  phases, 
and  on  the  acceptance  and  termination  criteria  will  be  explained  in 
detail  in  the  following  subsections. 

4.2.  A randomized  construction  procedure  for  the  generation 
of  initial  solutions 

Each  initial  solution  is  generated  by  means  of  the  following  ran- 
domized construction  procedure.  At  the  beginning,  all  (positive)  sup- 
plies are  sequenced  randomly.  According  to  this  sequence,  one  sup- 
ply after  another  is  assigned  to  a vehicle.  Starting  from  the  vehicle 
with  the  lowest  index,  it  is  checked  whether  the  supply  can  (still)  be 
accommodated  by  the  respective  vehicle.  This  is  the  case  if  (i)  a com- 
partment has  already  been  opened  for  the  respective  product  type 
(because  a supply  of  this  product  type  has  already  been  assigned 
to  the  vehicle)  or  it  still  can  be  opened  for  this  product  type  (be- 
cause there  is  still  enough  vehicle  capacity  available  and  the  maxi- 
mum number  of  compartments  is  not  exceeded)  and  (ii)  the  size  of 
the  compartment  is  at  least  as  large  as  the  supply.  If  the  conditions  (i) 
and  (ii)  hold,  then  the  supply  is  assigned  to  the  vehicle;  otherwise  the 
next  vehicle  will  be  checked.  The  procedure  stops  when  all  supplies 
have  been  assigned.  We  remark  that  the  assignment  of  supplies  to 
vehicles  also  constitutes  an  assignment  of  locations  to  vehicles.  The 
latter  is  used  as  a basis  for  the  determination  of  an  initial  set  of  routes 
by  application  of  the  Lin-Kernighan  Heuristic  (Lin  & Kernighan, 
1973). 

We  further  note  that  the  described  procedure  of  assigning  sup- 
plies to  vehicles  may  result  in  an  assignment,  which  contains  more 
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input:  problem  data,  number  of  neighborhood  structures  kmax; 

-,best 

Z : = oo ; 

do 

generate  an  initial  solution  x with  objective  function  value  Z(x)  randomly; 
set  neighborhood  structure  index  k to  k:=  1; 

do 

select  a neighbor  xJ  from  neighborhood  structure  k of  x randomly  (shaking); 
apply  local  search  to  x3  and  determine  a solution  x3  3 (improvement  phase  #1); 
if  Z(x3 3 ) satisfies  the  acceptance  criterion  then 
x : = x3  3 ; k:=  1; 
if  Z(x)  < Zbest  then 

best 

x : = x; 

apply  local  search  procedure  to  xbest  (improvement  phase  #2); 

Ztest:=  Z(xbest) ; 

endif 

else 

k:=  k + 1; 
if  k = kmax  + 1 then 
k:=  1; 
endif 
endif 

until  vns  termination  criterion  is  satisfied 
until  ms  termination  criterion  is  satisfied 
output:  xbest; 


Fig.  2.  Pseudo-code  of  a VNS  for  the  MCVRP-FCS. 


vehicles  than  actually  available.  In  this  case,  a large  penalty  is  added 
to  the  objective  function  for  each  additional  vehicle.  Consequently, 
the  algorithm  is  directed  to  finding  solutions  with  a feasible  num- 
ber of  vehicles  during  the  improvement  phases.  Once  a solution  with 
a feasible  number  of  vehicles  has  been  found,  then  only  solutions 
feasible  with  respect  to  the  number  of  available  vehicles  will  be 
generated. 

4.3.  Solution  space 

Throughout  the  search  of  the  solution  space,  not  only  feasible  so- 
lutions but  also  solutions  which  are  infeasible  with  regard  to  the  ve- 
hicle capacity  can  be  accepted  as  incumbent  solutions.  In  order  to 
deal  with  these  infeasibilities,  a penalty  term  is  added  to  the  total 
cost  of  a solution  within  the  objective  function  (Gendreau,  Hertz  and 
Laporte,  1994).  This  penalty  term  is  dependent  on  the  extent  accord- 
ing to  which  the  vehicle  capacities  are  violated  (excess  capacity  uti- 
lization) and  a penalty  factor  a.  For  the  determination  of  the  excess 
capacity  utilization  Q£  of  a vehicle  k,  the  selected  compartment  sizes 
are  considered  instead  of  the  actual  supplies  transported  in  these 
compartments: 

Q/[  = max  Qmy  pkm  Q. 

y peP  meM 

The  modified  objective  function  value  Z0  is  then  obtained  by 
adding  the  weighted  sum  of  the  total  excess  capacity  utilization  of 
all  vehicles  to  the  original  objective  function  value  Z0: 

Z0  = Z0  Q i< 

keK 

The  value  of  the  penalty  factor  a is  adjusted  dynamically.  At  the 
beginning  of  each  loop,  the  factor  is  set  to  an  identical  initial  value.  If 
no  feasible  solution  was  determined  for  a certain  number  of  consec- 
utive iterations,  the  factor  is  increased.  If  a feasible  solution  has  even- 
tually been  found,  the  factor  is  reset  to  its  initial  value.  In  this  man- 
ner, infeasibilities  get  penalized  only  lightly  during  the  early  stages 
of  the  application  of  the  algorithm.  However,  the  longer  the  solutions 
remain  infeasible,  the  stronger  the  penalization  gets.  Consequently, 


the  algorithm  is  directed  to  move  to  a feasible  region  of  the  solution 
space. 

4.4.  Neighborhood  structures 

The  crucial  part  of  the  design  of  a VNS  approach  consists  in  defin- 
ing the  neighborhood  structures  from  which  the  solutions  are  to  be 
selected  randomly  in  the  shaking  step.  Ten  different  neighborhood 
structures  have  been  implemented  in  order  to  deal  with  the  vari- 
ous decisions,  which  have  to  be  made  when  solving  the  MCVRP-FCS. 
They  can  be  distinguished  into  supply-related,  location-related,  and 
product  type-related  neighborhood  structures,  focusing  on  changes 
of  supply-vehicle-assignments,  location-vehicle-assignments,  and 
product-vehicle-assignments,  respectively. 

Two  supply-related  neighborhood  structures  are  used,  of  which 
the  first  one  is  based  on  supply  shifts,  while  the  second  one  is  based 
on  supply  swaps.  A supply  shift  selects  a single  supply  randomly, 
deletes  it  from  its  current  vehicle  and  inserts  it  into  another,  ran- 
domly selected,  vehicle.  A supply  swap  selects  two  individual  sup- 
plies randomly  from  two  different  vehicles  and  exchanges  both  sup- 
plies. Through  application  of  these  and  all  following  neighborhood 
structures,  it  might  also  be  necessary  to  adjust  routes.  These  adjust- 
ments are  performed  in  the  following  ways:  (1 ) In  case  all  supplies  at 
a location  are  eliminated  from  one  vehicle,  this  location  is  deleted 
from  the  corresponding  route  and  the  edge  between  the  locations 
which  were  adjacent  to  the  eliminated  location  is  added;  (2)  in  case 
a location  has  to  be  added  to  a vehicle,  it  is  inserted  at  a random  po- 
sition in  the  respective  route.  For  both  neighborhood  structures  an 
example  is  shown  in  Fig.  3.  Both  examples  are  based  on  the  problem 
instance  shown  in  Section  2.  For  the  purpose  of  illustration,  three  ve- 
hicles are  used  in  the  initial  solution  instead  of  two  vehicles  as  in 
Fig.  1. 

Four  different  location-related  neighborhood  structures  have 
been  implemented,  based  on  complete  location  shifts,  partial  lo- 
cation shifts,  location  splits,  and  location  swaps.  For  a complete 
location  shift,  one  location-vehicle-assignment  from  the  incumbent 
solution  is  randomly  selected.  Then,  all  supplies  related  to  this  spe- 
cific location,  which  are  currently  assigned  to  the  vehicle,  are  moved 
collectively  to  another  vehicle.  In  a partial  location  shift,  only  a 
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a)  Supply  shift:  The  supply  of  product  type  R at  location  1 
is  moved  from  vehicle  1 to  vehicle  2. 


b)  Supply  swap:  The  supply  of  product  type  Y at  location  1 
and  the  supply  of  product  type  R at  location  2 are 
exchanged  between  vehicles  1 and  2. 


Fig.  3.  Illustration  of  supply-related  neighborhood  structures. 


randomly  determined  subset  of  these  supplies  is  shifted.  A location 
split  operator  works  in  a similar  way  as  a complete  location  shift  op- 
erator except  for  the  difference  that  each  supply  may  be  shifted  to  a 
different  vehicle.  In  a location  swap,  two  location-vehicle- 
assignments  are  selected  randomly  and  the  corresponding  supplies 
of  the  locations  are  exchanged  between  the  vehicles.  Pretests  have 
also  shown  that  it  is  sensible  to  change  more  than  one  location- 
vehicle-assignment  simultaneously,  e.g.  for  a complete  location  shift 
not  only  one  location  but  two  locations  are  shifted  to  other  vehicles. 
The  four  described  location-related  neighborhood  structures  are 
illustrated  in  Fig.  4. 

Four  different  options  for  the  product  type-related  neighborhood 
structures  have  been  implemented.  Analogously  to  the  location- 
related  neighborhood  structures,  they  are  established  by  complete 
product  shifts,  partial  product  shifts,  product  splits,  and  product 
swaps.  In  contrast  to  location-related  operators,  product-vehicle- 
assignments  instead  of  location-vehicle-assignments  are  selected 
randomly.  As  for  a product  shift,  the  supplies  of  the  selected  product 
type  are  moved  either  completely  (complete  shift)  or  partially  (par- 
tial shift)  to  another  vehicle.  For  the  product  split,  the  supplies  may 
be  assigned  to  different  vehicles.  In  a product  swap,  two  product- 
vehicle-assignments  are  selected  randomly  and  the  corresponding 
total  supplies  are  exchanged  between  the  two  vehicles.  An  illus- 
tration of  the  product-related  neighborhood  structures  is  shown  in 
Fig.  5. 

In  the  VNS  algorithm  for  the  MCVRP-FCS,  the  described  neighbor- 
hood structures  have  been  implemented  in  the  following  sequence: 
demand  shift,  demand  swap,  partial  location  shift,  complete  loca- 
tion shift,  location  swap,  location  split,  partial  product  shift,  complete 
product  shift,  product  swap,  and  product  split.  Only  such  moves  are 
considered  which  neither  violate  the  maximum  number  of  product 
types  per  vehicle  nor  the  number  of  available  vehicles. 

4.5.  Improvement  phases 

During  the  improvement  phase  immediately  following  the  shak- 
ing step  (improvement  phase  #1 ),  all  routes  of  a solution  which  have 
been  changed  by  the  shaking  step  are  improved  by  a local  search 
procedure.  This  procedure  uses  the  3-edge  neighborhood  structure 


proposed  by  Lin  (1965),  where  the  neighborhood  is  established  by 
an  operator,  according  to  which  three  edges  are  eliminated  from 
a route  and  three  new  edges  are  inserted.  The  search  applies  the 
first  improvement  principle,  i.e.  the  first  route  is  accepted  from  a 
neighborhood  which  has  smaller  route  costs  than  the  current  one. 
The  improvement  phase  terminates  when  no  further  improved  route 
can  be  identified  in  the  neighborhood  of  the  currently  considered 
route. 

In  the  improvement  phase,  which  follows  the  identification  of  a 
new  best  solution  (improvement  phase  #2),  the  well-known  Lin- 
Kernighan  Fleuristic  (Lin  & Kernighan,  1973)  is  applied  to  each  of  the 
routes  in  the  solution.  For  our  implementation  we  used  the  code  pro- 
vided by  Flelsgaun  (2000). 

4.6.  Acceptance  criterion 

VNS  algorithms  tend  to  get  stuck  in  local  optima  (cf.  Flansen  & 
Mladenovic,  2001).  In  order  to  avoid  this  drawback,  we  do  not  only 
use  several  starts  of  the  heuristic  but  also  an  acceptance  criterion 
which  allows  for  accepting  solutions  with  objective  function  values 
worse  than  those  of  the  best  solution  found  in  a loop.  The  accep- 
tance criterion  used  for  this  algorithm  is  based  on  threshold  accept- 
ing (Dueck  & Scheuer,  1990)  with  a dynamically  self-adjusting  thresh- 
old factor. 

The  dynamic  self-adjustment  principle  is  similar  to  the  adjust- 
ment procedure  for  the  penalty  factor  described  in  Section  4.3.  The 
factor  starts  with  an  initial  value  of  1.  After  a certain  number  of 
consecutive  iterations  in  which  no  new  incumbent  solution  was  ac- 
cepted, this  factor  is  increased  by  a certain  value.  If  a solution  is 
eventually  accepted  as  incumbent  solution,  the  factor  is  reset  to  1 
again.  In  this  way,  the  threshold  increases  slowly  if  the  algorithm  is 
not  able  to  find  a better  solution  in  the  region  of  the  solution  space 
which  is  currently  explored,  enabling  the  algorithm  to  overcome  local 
optima. 

4.7.  Termination  criteria 

A loop  terminates  after  a certain  number  of  iterations  without 
finding  a new  best  solution.  The  algorithm  terminates  after  a certain 
number  of  loops  without  finding  a new  globally  best  solution. 
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a)  Complete  location  shift:  All  supplies  at  location  1 which 
are  currently  assigned  to  vehicle  1 are  moved  collectively  to 
vehicle  3 . 
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b)  Partial  location  shift:  Some  supplies  at  location  1 which 
are  currently  assigned  to  vehicle  1 are  moved  collectively  to 
vehicle  3 . 


c)  Location  swap:  All  supplies  at  location  1 which  are 
currently  assigned  to  vehicle  1 and  all  supplies  at  location  3 
which  are  currently  assigned  to  vehicle  3 are  exchanged 
between  vehicles  1 and  3 . 


d)  Location  split:  All  supplies  at  location  1 which  are 
currently  assigned  to  vehicle  1 are  moved  individually  to 
vehicles  2 and  3 . 


Fig.  4.  Illustration  of  location-related  neighborhood  structures. 


5.  Numerical  experiments 

5 A.  Overview 

The  above-presented  model  of  the  MCVRP-FCS  was  imple- 
mented - with  a separation  procedure  for  the  subtour-elimination 
constraints  - in  Microsoft  Visual  C++  2012  using  CPLEX  12.5.  For  the 
MS-VNS  algorithm,  Microsoft  Visual  C++  2012  has  been  used. 

Three  sets  of  numerical  experiments  were  performed  in  order 
to  evaluate  the  two  approaches.  In  the  first  set,  the  exact  solution 
approach  and  the  VNS  algorithm  have  been  applied  to  relatively 
small,  randomly-generated  problem  instances.  These  experiments 
serve  two  purposes.  Firstly,  they  provide  an  impression  of  the  limits 
of  the  problem  size  up  to  which  the  MCVRP-FCS  can  be  solved  opti- 
mally by  means  of  a standard  LP  solver  in  reasonable  computing  time. 
Secondly,  by  comparing  the  results  from  the  MS-VNS  approach  to  the 


results  obtained  by  the  exact  approach,  the  solution  quality  of  the 
MS-VNS  approach  can  be  assessed. 

In  a second  set  of  experiments,  large,  randomly-generated  in- 
stances of  the  MCVRP-FCS  were  considered,  to  which  the  MS-VNS 
algorithm  was  applied.  These  experiments  have  been  performed  in 
order  to  study  the  behavior  of  the  algorithm  in  greater  detail.  They 
were  meant  to  provide  insights  into  how  different  problem  param- 
eters affect  solution  quality  and  computing  times  and,  in  particular, 
which  parameters  make  the  problem  difficult  to  solve. 

In  the  third  set  of  experiments,  the  MS-VNS  algorithm  was  applied 
to  problem  instances  from  practice.  These  experiments  were  meant 
to  assess  the  benefits  from  introducing  vehicles  with  compartments 
of  flexible  size  in  comparison  to  the  utilization  of  vehicles  with  a sin- 
gle compartment  only. 

All  experiments  were  performed  on  a 3.2  gigahertz  and  8 giga- 
bytes RAM  personal  computer. 
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a)  Complete  product  shift:  All  supplies  of  product  type  Y 
which  are  currently  assigned  to  vehicle  1 are  moved 
collectively  to  vehicle  2. 


b)  Partial  product  shift:  Some  supplies  of  product  type  Y 
which  are  currently  assigned  to  vehicle  1 are  moved 
collectively  to  vehicle  2. 
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c)  Product  swap:  All  supplies  of  product  type  Y which  are 
currently  assigned  to  vehicle  1 and  all  supplies  of  product 
type  B which  are  currently  assigned  to  vehicle  2 are 
exchanged  between  vehicles  1 and  2. 


d)  Product  split:  All  supplies  of  product  type  Y which  are 
currently  assigned  to  vehicle  1 are  moved  individually  to 
vehicles  2 and  3 . 


Fig.  5.  Illustration  of  product-related  neighborhood  structures. 


5.2.  Problem  generation 

As  we  have  explained  above,  the  MCVRP-FCS  is  different  from 
other  multi-compartment  vehicle  routing  problems  described  in  the 
literature  so  far.  Therefore,  we  could  not  use  benchmark  instances 
from  the  literature  in  our  experiments.  Instead,  we  randomly  gen- 
erated new  problem  instances. 

In  order  to  establish  different  classes  of  problem  instances,  the  pa- 
rameters of  the  MCVRP-FCS  have  been  dealt  with  in  the  following 
way.  The  number  of  locations  n,  the  number  of  product  types  |P|,  and 
the  number  of  compartments  m have  been  taken  as  controllable  pa- 
rameters. The  respective  values  chosen  for  the  experiments  are  given 
below.  Concerning  the  number  of  product  types  of  which  supplies 
are  available  at  each  location  (number  of  supplies),  three  options  for 
a supply  parameter  s were  considered:  small,  medium,  and  large.  A 
small  number  of  supplies  corresponds  to  a situation  in  which  each 
location  provides  supplies  for  at  least  one  product  type  and  at  most 


for  one  third  of  the  product  types  (s  = 1,  small  number  of  supplies). 
A medium  number  of  supplies  represents  a situation  where  supplies 
of  at  least  more  than  one  third  and  at  most  two  thirds  of  the  product 
types  are  available  at  each  location  (s  = 2,  medium  number  of  sup- 
plies). Finally,  for  a large  number  of  supplies,  supplies  are  available  at 
each  location  of  at  least  more  than  two  thirds  of  the  product  types  up 
to  \P\  (i.e.  all)  product  types  (s  = 3,  large  number  of  supplies). 

The  actual  supplies  for  all  product  types  at  all  locations  sip  were 
generated  by  means  of  the  following  procedure.  At  first,  an  initial 
total  supply  (STotal)  for  all  product  types  across  all  locations  was 
randomly  generated  from  the  integer  interval  [125  • n;  175  • n].  This 
total  supply  was  then  split  into  individual  supplies  at  the  locations  in 
several  steps.  Firstly,  based  on  the  respective  value  of  the  parameter 
s,  a number  for  the  supplied  product  types  (si,  i e V\{0})  was  gener- 
ated randomly  for  each  location  i from  the  corresponding  parameter 
interval  for  s.  In  order  to  determine  which  specific  product  types  at  a 
specific  location  have  an  actual  supply,  i.e.  sip  > 0,  a random  number 
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out  of  the  interval  [0,  0.5]  for  the  first  third  of  product  types,  [0,  0.3] 
for  the  second  third  of  product  types,  and  [0,  0.2]  for  the  last  third  of 
product  types  was  drawn.  For  each  location  z,  these  random  numbers 
were  sorted  in  decreasing  order,  whereas  the  first  st  product  types 
in  this  sequence  were  selected  to  have  a supply  greater  than  zero. 
Then,  the  total  supply  quantities  for  each  product  type  (Sjotal,  p e P) 
were  determined  in  such  a manner  that  50  percent  of  the  total  supply 
STotal  was  assigned  to  the  first  third  of  the  product  types,  30  percent 
of  the  total  supply  to  the  second  third,  and  20  percent  of  the  total 
supply  to  the  last  third  of  the  product  types.  Finally,  the  actual  sup- 
plies (SjP,  i e V\{0},  p e P,  sip  > 0)  were  generated  by  drawing  a ran- 
dom number  for  each  supply,  which  were  subsequently  normalized 
according  to  the  total  supply  of  the  corresponding  product  type  S]°tal. 
The  biased  distribution  of  product  type  supplies  was  used  in  order  to 
generate  realistic  instances,  since  in  the  glass  collection  problem  the 
different  glass  types  are  more  and  less  frequently  disposed. 

Finally,  the  vehicle  capacity  was  set  to  1000  for  each  instance  and 
the  number  of  vehicles  was  determined  by  solving  a corresponding 
bin-packing  problem  exactly  in  order  to  guarantee  a feasible  solution 
to  the  problem  instance. 

All  instances  can  be  found  online  at  http://www.mansci.ovgu. 
de/mansci/en/Research/Materials/2014-p-394.html. 

5.3.  Experiments  with  small , randomly-generated  instances 

The  following  parameter  values  were  used  for  the  generation  of 
instances  for  the  first  set  of  experiments: 

• number  of  customer  locations  n = 10; 

• number  of  products  types  \P\  e {3,  6,  9}; 

• number  of  compartments  m e {2,  3}  for  \P\  = 3,  m e {2,  4,  6} 

for  |P|  = 6,  me  {2, 4,  7,  9}  for  |P|  = 9; 

• supply  parameter  s e {1,  2,  3}. 


that  the  MCVR-FCS  becomes  more  difficult  to  solve  (i.e.  the  re- 
spective computing  times  increase)  with  an  increasing  number 
of  product  types,  with  an  increasing  number  of  supplies  and 
with  a decreasing  number  of  compartments.  Especially  the  lat- 
ter observation  demonstrates  that  the  specific  (partial)  decision 
of  this  problem,  namely  the  assignment  of  product  types  to  ve- 
hicles, has  a significant  impact  on  the  complexity  of  the  prob- 
lem and  makes  it  different  from  other  vehicle  routing  problems. 
We  further  note  that  - with  respect  to  the  computing  times  - 
the  exact  approach  presented  here  cannot  be  expected  to  be  com- 
petitive for  large,  real-world  instances. 

Regarding  the  solution  quality,  it  can  be  observed  that  the 
MS-VNS  algorithm  performs  very  well  on  the  instances  from  these 
problem  classes.  From  the  1350  instances,  the  MS-VNS  algorithm 
solves  1334  (98.8  percent)  optimally.  The  average  deviation  from  the 
optimal  objective  function  value  (across  all  instances)  only  amounted 
to  0.02  percent.  This  is  particularly  remarkable  since  the  average 
computing  time  per  problem  instance  needed  by  the  MS-VNS  algo- 
rithm (12.97  seconds)  represents  only  3.0  percent  of  the  average  com- 
puting time  needed  by  the  exact  approach  (431.31  seconds).  Regard- 
ing the  maximal  computing  times,  this  observation  is  even  more  con- 
siderable, since  the  maximal  computing  time  for  the  exact  approach 
is  21  hours  whereas  the  maximal  computing  time  for  the  MS-VNS  is 
only  53  seconds. 

As  for  the  computing  times  of  the  MS-VNS  algorithm,  it  can  be 
observed  that  they  increase  with  the  number  of  product  types  and 
the  number  of  demands.  In  contrast  to  the  exact  approach,  the  impact 
of  the  maximum  number  of  compartments  on  the  computing  time  is 
inconsistent.  For  instances  with  three  or  six  product  types,  no  pattern 
is  observable;  for  instances  with  nine  product  types,  the  computing 
times  decrease  with  increasing  number  of  available  compartments. 

5.4.  Experiments  with  large , randomly-generated  instances 


The  number  of  locations  considered  in  this  set  of  experiments  may 
appear  to  be  small  at  first  sight,  but  it  has  to  be  emphasized  that, 
apart  from  the  number  of  locations,  also  the  total  number  of  supplies 
determines  the  problem  size.  Since  the  number  of  different  supplies 
at  each  location  and  the  maximum  number  of  compartments  per  ve- 
hicle make  the  MCVRP-FCS  special,  we  have  chosen  to  explore  these 
parameters  in  the  first  place  at  the  expense  of  considering  not  more 
than  10  customer  locations.  This  confinement  was  necessary  in  order 
to  limit  the  respective  optimization  models  to  sizes  which  were  still 
manageable  by  CPLEX. 

The  parameter  settings  resulted  in  27  problem  classes,  for  each 
of  which  50  instances  have  been  generated.  To  all  1350  instances, 
the  exact  approach  and  the  MS-VNS  algorithm  have  been  applied. 
500n|P|  iterations  without  improvement  (vns  termination  criterion; 
cf.  Fig.  2)  and  (n|P|)/6  loops  without  improvement  (ms  termination 
criterion)  were  used  as  termination  criteria  for  the  MS-VNS  algo- 
rithm. For  each  problem  class,  Table  1 presents  the  number  of  in- 
stances, which  have  been  solved  to  a proven  optimum  by  the  ex- 
act approach  (#opt.exact)  and  by  the  MS-VNS  algorithm  (#opt.vns), 
the  average  objective  function  value  (total  cost)  per  problem  instance 
obtained  by  the  exact  approach  (tc.exact)  and  by  the  MS-VNS  algo- 
rithm (tc.vns),  and  computing  times  per  problem  instance  in  sec- 
onds for  the  exact  approach  (cpu.exact)  and  the  MS-VNS  algorithm 
(cpu.vns).  For  the  exact  approach,  also  the  maximal  computing  time 
(cpu.exactmax)  needed  for  an  instance  is  reported  for  each  problem 
class.  Furthermore,  the  average  per-instance  deviation  (aver.gap)  of 
the  objective  function  value  obtained  by  the  MS-VNS  algorithm  from 
the  objective  function  value  obtained  by  the  exact  approach,  and  the 
corresponding  maximum  deviation  (max.gap)  across  all  instances  in 
a class  are  presented. 

The  exact  approach  managed  to  solve  all  of  the  1350  in- 
stances with  a maximal  computing  time  of  21  hours.  It  appears 


For  the  second  set  of  experiments,  a set  of  large  instances  was 
generated.  These  instances  consider  50  locations  and  up  to  9 product 
types,  which  results  in  a maximum  of  450  demands.  For  the  number 
of  product  types,  the  number  of  compartments,  and  the  number  of 
demands,  the  same  combinations  of  parameters  as  in  the  first  set  of 
experiments  were  used.  For  each  combination,  one  instance  was  ran- 
domly generated.  Based  on  these  experiments,  the  solution  quality  of 
the  MS-VNS  was  evaluated. 

Since  we  were  neither  able  to  use  benchmarks  from  the  litera- 
ture nor  to  solve  any  large  instance  to  optimality,  we  adopted  the  fol- 
lowing approach.  Each  instance  was  solved  for  360  minutes  by  the 
MS-VNS.  The  best  solutions  obtained  after  10,  20,  30,  40,  50,  and  60, 
120, 180, 240, 300,  and  360  minutes  were  recorded  and  the  respective 
objective  function  values  were  compared  to  the  best  solution  found 
after  360  minutes.  Again,  500n|P|  iterations  without  improvement 
per  loop  were  used. 

Table  2 shows  for  each  instance  the  total  costs  of  the  solutions 
found  (tc)  during  the  first  hour  of  computation  in  steps  of  10  minutes 
and  the  corresponding  deviations  to  the  best  known  solution  (gap 
w.r.t.  BKS).  Moreover,  Table  3 lists  similar  results  for  the  6 hours  of 
computation  in  steps  of  60  minutes.  For  each  instance,  the  objective 
function  value  of  the  best  known  solution  is  indicated  by  BKS  in  both 
tables. 

Table  2 shows  that,  on  average,  solutions  found  after  10  min- 
utes have  an  objective  function  value  which  deviates  5.62  percent 
from  the  best  solution  found  after  360  minutes.  After  60  minutes 
of  running  time,  this  gap  is  decreased  to  2.03  percent  on  average. 
Furthermore,  one  can  observe  a decrease  in  the  change  of  the  gaps: 
Whereas  the  solutions  were  improved  by  1.63  percent  and  1.04  per- 
cent on  average  during  the  second  and  third  10-minute  intervals,  re- 
spectively, this  change  decreased  to  0.19  percent,  0.27  percent,  and 
0.45  percent  for  the  remaining  three  10-minute  intervals.  Similar 
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Table  1 

Results  for  small  instances. 


Parameters 

cplex 

vns 

gap  w.r.t.  opt 

\P\ 

m 

s 

#opt.exact 

tc.exact 

cpu.exact 

cpu.exactmax 

#opt.vns 

tc.vns 

cpu.vns 

avg.gap  (percent) 

max.gap  (percent) 

3 

2 

1 (small) 

50 

366.9 

5.43 

52.61 

50 

366.9 

1.19 

0.00 

0.00 

2 (medium) 

50 

465.3 

11.43 

194.87 

50 

465.35 

1.54 

0.00 

0.00 

3 (large) 

50 

596.2 

101.33 

4686.80 

50 

596.16 

2.02 

0.00 

0.00 

3 

1 (small) 

50 

349.1 

2.60 

25.84 

50 

349.08 

1.28 

0.00 

0.00 

2 (medium) 

50 

337.1 

2.66 

33.65 

50 

337.15 

1.64 

0.00 

0.00 

3 (large) 

50 

336.6 

5.00 

94.45 

50 

336.64 

1.89 

0.00 

0.00 

6 

2 

1 (small) 

50 

550.0 

5.91 

28.53 

50 

549.98 

3.80 

0.00 

0.00 

2 (medium) 

50 

767.5 

4.83 

21.26 

50 

767.48 

5.47 

0.00 

0.00 

3 (large) 

50 

919.8 

1.81 

3.96 

50 

919.78 

6.44 

0.00 

0.00 

4 

1 (small) 

50 

367.0 

5.04 

47.89 

50 

366.98 

5.58 

0.00 

0.00 

2 (medium) 

50 

492.1 

27.79 

120.60 

49 

492.13 

7.60 

0.00 

0.21 

3 (large) 

50 

583.8 

326.15 

2895.40 

49 

584.46 

9.94 

0.15 

7.39 

6 

1 (small) 

50 

352.0 

4.05 

72.42 

50 

351.98 

6.03 

0.00 

0.00 

2 (medium) 

50 

358.0 

5.94 

27.14 

50 

357.98 

7.73 

0.00 

0.00 

3 (large) 

50 

357.4 

8.96 

62.29 

50 

357.42 

8.44 

0.00 

0.00 

9 

2 

1 (small) 

50 

724.2 

502.44 

13866.00 

49 

724.37 

18.00 

0.03 

1.39 

2 (medium) 

50 

1114.3 

1921.59 

12333.90 

50 

1114.25 

25.68 

0.00 

0.00 

3 (large) 

50 

1419.1 

5400.85 

76049.40 

50 

1419.09 

36.19 

0.00 

0.00 

4 

1 (small) 

50 

463.0 

12.87 

76.28 

50 

463.05 

15.67 

0.00 

0.00 

2 (medium) 

50 

647.6 

267.69 

4653.76 

50 

647.57 

22.87 

0.00 

0.00 

3 (large) 

50 

815.3 

1324.58 

4712.93 

50 

815.27 

26.46 

0.00 

0.00 

7 

1 (small) 

50 

360.2 

5.25 

60.90 

50 

360.16 

17.47 

0.00 

0.00 

2 (medium) 

50 

456.3 

93.21 

2239.98 

44 

457.40 

22.15 

0.22 

6.95 

3 (large) 

50 

564.5 

1550.47 

10854.80 

46 

565.32 

30.00 

0.15 

3.22 

9 

1 (small) 

50 

355.0 

3.94 

18.95 

50 

354.98 

17.65 

0.00 

0.00 

2 (medium) 

50 

392.5 

10.53 

74.96 

50 

392.48 

22.54 

0.00 

0.00 

3 (large) 

50 

383.4 

32.94 

459.02 

47 

383.81 

25.01 

0.12 

4.04 

Minimum 

50 

336.6 

1.81 

3.96 

44 

336.64 

1.19 

0.00 

0.00 

Average 

50 

551.6 

431.31 

4954.39 

49.4 

551.75 

12.97 

0.02 

0.86 

Maximum 

50 

1419.1 

5400.85 

76049.40 

50 

1419.09 

36.19 

0.22 

7.39 
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Table  2 

Results  for  large  instances,  1 hour  of  computing  time. 
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\p\ 

in 

s 

Best  solution  found  after 

BKS 

10  minutes 

20  minutes 

30  minutes 

40  minutes 

50  minutes 

60  minutes 

tc 

gap  w.r.t.  BKS  tc 
(percent) 

gap  w.r.t.  BKS  tc 
(percent) 

gap  w.r.t.  BKS  tc 
(percent) 

gap  w.r.t.  BKS  tc 
(percent) 

gap  w.r.t.  BKS  tc 
(percent) 

gap  w.r.t.  BKS 
(percent) 

3 

2 

1 (small) 

1117.79 

3.96 

1117.79 

3.96 

1101.82 

2.47 

1097.24 

2.05 

1097.24 

2.05 

1080.68 

0.51 

1075.21 

2 (medium) 

1136.31 

0.00 

1136.31 

0.00 

1136.31 

0.00 

1136.31 

0.00 

1136.31 

0.00 

1136.31 

0.00 

1136.31 

3 (large) 

1529.73 

4.84 

1459.09 

0.00 

1459.09 

0.00 

1459.09 

0.00 

1459.09 

0.00 

1459.09 

0.00 

1459.09 

3 

1 (small) 

1347.26 

6.96 

1259.55 

0.00 

1259.55 

0.00 

1259.55 

0.00 

1259.55 

0.00 

1259.55 

0.00 

1259.55 

2 (medium) 

1792.21 

10.89 

1792.21 

10.89 

1667.96 

3.20 

1667.96 

3.20 

1667.96 

3.20 

1667.96 

3.20 

1616.20 

3 (large) 

1830.68 

2.75 

1830.68 

2.75 

1830.68 

2.75 

1830.68 

2.75 

1830.68 

2.75 

1830.68 

2.75 

1781.74 

6 

2 

1 (small) 

1323.69 

0.00 

1323.69 

0.00 

1323.69 

0.00 

1323.69 

0.00 

1323.69 

0.00 

1323.69 

0.00 

1323.69 

2 (medium) 

1970.85 

8.07 

1970.85 

8.07 

1844.73 

1.15 

1844.73 

1.15 

1844.73 

1.15 

1844.73 

1.15 

1823.69 

3 (large) 

2879.94 

20.16 

2396.85 

0.00 

2396.85 

0.00 

2396.85 

0.00 

2396.85 

0.00 

2396.85 

0.00 

2396.85 

4 

1 (small) 

1112.26 

0.00 

1112.26 

0.00 

1112.26 

0.00 

1112.26 

0.00 

1112.26 

0.00 

1112.26 

0.00 

1112.26 

2 (medium) 

1627.56 

6.42 

1627.56 

6.42 

1627.56 

6.42 

1627.56 

6.42 

1627.56 

6.42 

1589.01 

3.90 

1529.42 

3 (large) 

2486.69 

8.06 

2469.11 

7.30 

2469.11 

7.30 

2434.65 

5.80 

2372.35 

3.09 

2372.35 

3.09 

2301.20 

6 

1 (small) 

1648.06 

3.48 

1647.35 

3.43 

1631.44 

2.44 

1631.44 

2.44 

1631.44 

2.44 

1631.44 

2.44 

1592.65 

2 (medium) 

2069.92 

3.48 

2069.92 

3.48 

2069.92 

3.48 

2069.92 

3.48 

2069.92 

3.48 

2052.98 

2.63 

2000.30 

3 (large) 

2472.03 

6.59 

2472.03 

6.59 

2472.03 

6.59 

2428.83 

4.73 

2387.18 

2.93 

2387.18 

2.93 

2319.24 

9 

2 

1 (small) 

1607.41 

0.71 

1607.41 

0.71 

1597.10 

0.06 

1597.10 

0.06 

1597.10 

0.06 

1597.10 

0.06 

1596.12 

2 (medium) 

2726.08 

2.95 

2702.56 

2.06 

2702.56 

2.06 

2702.56 

2.06 

2702.56 

2.06 

2702.56 

2.06 

2648.00 

3 (large) 

3343.61 

7.10 

3268.01 

4.68 

3268.01 

4.68 

3268.01 

4.68 

3268.01 

4.68 

3268.01 

4.68 

3121.83 

4 

1 (small) 

1480.58 

8.33 

1480.58 

8.33 

1463.77 

7.10 

1463.77 

7.10 

1463.47 

7.07 

1393.81 

1.98 

1366.79 

2 (medium) 

1541.96 

0.00 

1541.96 

0.00 

1541.96 

0.00 

1541.96 

0.00 

1541.96 

0.00 

1541.96 

0.00 

1541.96 

3 (large) 

1977.32 

0.00 

1977.32 

0.00 

1977.32 

0.00 

1977.32 

0.00 

1977.32 

0.00 

1977.32 

0.00 

1977.32 

7 

1 (small) 

1601.54 

13.64 

1534.57 

8.89 

1509.99 

7.14 

1490.15 

5.73 

1490.15 

5.73 

1490.15 

5.73 

1409.34 

2 (medium) 

2517.56 

7.32 

2517.56 

7.32 

2494.68 

6.34 

2494.68 

6.34 

2494.68 

6.34 

2494.68 

6.34 

2345.86 

3 (large) 

3830.43 

6.56 

3711.65 

3.25 

3711.65 

3.25 

3711.65 

3.25 

3612.92 

0.51 

3612.92 

0.51 

3594.65 

9 

1 (small) 

1843.51 

6.36 

1843.51 

6.36 

1789.81 

3.26 

1789.81 

3.26 

1789.81 

3.26 

1789.81 

3.26 

1733.29 

2 (medium) 

2816.5 

7.21 

2816.50 

7.21 

2793.74 

6.35 

2793.74 

6.35 

2793.74 

6.35 

2793.74 

6.35 

2627.05 

3 (large) 

4117.89 

5.83 

4117.89 

5.83 

4023.72 

3.41 

4023.72 

3.41 

4023.72 

3.41 

3942.29 

1.32 

3890.97 

Minimum 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Average 

5.62 

3.98 

2.94 

2.75 

2.48 

2.03 

Maximum 

20.16 

10.89 

7.30 

7.10 

7.07 

6.35 
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Table  3 

Results  for  large  instances,  6 hours  of  computing  time. 
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Best  solution  found  after 

60  minutes 

120  minutes 

180  minutes 

240  minutes 

300  minutes 

360  minutes 

tc 

gap  w.r.t.  BKS  tc 

(percent) 

gap  w.r.t.  BKS  tc 

(percent) 

gap  w.r.t.  y tc 

(percent) 

gap  w.r.t.  BKS  tc 

(percent) 

gap  w.r.t.  y tc  (BKS) 

(percent) 

gap  w.r.t.  BKS 
(percent) 

3 

2 

1 (small) 

1080.68 

0.51 

1080.68 

0.51 

1075.76 

0.05 

1075.21 

0.00 

1075.21 

0.00 

1075.21 

0.00 

2 (medium) 

1136.31 

0.00 

1136.31 

0.00 

1136.31 

0.00 

1136.31 

0.00 

1136.31 

0.00 

1136.31 

0.00 

3 (large) 

1459.09 

0.00 

1459.09 

0.00 

1459.09 

0.00 

1459.09 

0.00 

1459.09 

0.00 

1459.09 

0.00 

3 

1 (small) 

1259.55 

0.00 

1259.55 

0.00 

1259.55 

0.00 

1259.55 

0.00 

1259.55 

0.00 

1259.55 

0.00 

2 (medium) 

1667.96 

3.20 

1616.20 

0.00 

1616.20 

0.00 

1616.20 

0.00 

1616.20 

0.00 

1616.20 

0.00 

3 (large) 

1830.68 

2.75 

1781.74 

0.00 

1781.74 

0.00 

1781.74 

0.00 

1781.74 

0.00 

1781.74 

0.00 

6 

2 

1 (small) 

1323.69 

0.00 

1323.69 

0.00 

1323.69 

0.00 

1323.69 

0.00 

1323.69 

0.00 

1323.69 

0.00 

2 (medium) 

1844.73 

1.15 

1844.73 

1.15 

1844.73 

1.15 

1844.73 

1.15 

1844.73 

1.15 

1823.69 

0.00 

3 (large) 

2396.85 

0.00 

2396.85 

0.00 

2396.85 

0.00 

2396.85 

0.00 

2396.85 

0.00 

2396.85 

0.00 

4 

1 (small) 

1112.26 

0.00 

1112.26 

0.00 

1112.26 

0.00 

1112.26 

0.00 

1112.26 

0.00 

1112.26 

0.00 

2 (medium) 

1589.01 

3.90 

1589.01 

3.90 

1571.95 

2.78 

1571.95 

2.78 

1529.42 

0.00 

1529.42 

0.00 

3 (large) 

2372.35 

3.09 

2354.77 

2.33 

2301.20 

0.00 

2301.20 

0.00 

2301.20 

0.00 

2301.20 

0.00 

6 

1 (small) 

1631.44 

2.44 

1609.81 

1.08 

1592.65 

0.00 

1592.65 

0.00 

1592.65 

0.00 

1592.65 

0.00 

2 (medium) 

2052.98 

2.63 

2052.98 

2.63 

2039.76 

1.97 

2012.06 

0.59 

2000.30 

0.00 

2000.30 

0.00 

3 (large) 

2387.18 

2.93 

2368.14 

2.11 

2368.14 

2.11 

2368.14 

2.11 

2319.24 

0.00 

2319.24 

0.00 

9 

2 

1 (small) 

1597.100 

0.06 

1596.12 

0.00 

1596.12 

0.00 

1596.12 

0.00 

1596.12 

0.00 

1596.12 

0.00 

2 (medium) 

2702.56 

2.06 

2675.37 

1.03 

2675.37 

1.03 

2661.55 

0.51 

2648.00 

0.00 

2648.00 

0.00 

3 (large) 

3268.01 

4.68 

3268.01 

4.68 

3248.53 

4.06 

3121.83 

0.00 

3121.83 

0.00 

3121.83 

0.00 

4 

1 (small) 

1393.81 

1.98 

1393.81 

1.98 

1393.81 

1.98 

1393.81 

1.98 

1393.81 

1.98 

1366.79 

0.00 

2 (medium) 

1541.96 

0.00 

1541.96 

0.00 

1541.96 

0.00 

1541.96 

0.00 

1541.96 

0.00 

1541.96 

0.00 

3 (large) 

1977.32 

0.00 

1977.32 

0.00 

1977.32 

0.00 

1977.32 

0.00 

1977.32 

0.00 

1977.32 

0.00 

7 

1 (small) 

1490.15 

5.73 

1472.07 

4.45 

1409.34 

0.00 

1409.34 

0.00 

1409.34 

0.00 

1409.34 

0.00 

2 (medium) 

2494.68 

6.34 

2464.19 

5.04 

2370.69 

1.06 

2370.69 

1.06 

2345.86 

0.00 

2345.86 

0.00 

3 (large) 

3612.92 

0.51 

3612.92 

0.51 

3612.92 

0.51 

3612.92 

0.51 

3594.65 

0.00 

3594.65 

0.00 

9 

1 (small) 

1789.81 

3.26 

1765.17 

1.84 

1759.60 

1.52 

1741.41 

0.47 

1741.41 

0.47 

1733.29 

0.00 

2 (medium) 

2793.74 

6.35 

2709.40 

3.13 

2709.40 

3.13 

2709.40 

3.13 

2705.09 

2.97 

2627.05 

0.00 

3 (large) 

3942.29 

1.32 

3942.29 

1.32 

3890.97 

0.00 

3890.97 

0.00 

3890.97 

0.00 

3890.97 

0.00 

Minimum 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Average 

2.03 

1.40 

0.79 

0.53 

0.24 

0.00 

Maximum 

6.35 

5.04 

4.06 

3.13 

2.97 

0.00 
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Fig.  6.  Average  gap  with  respect  to  the  best  known  solution  plotted  against  computing 
time. 

results  are  shown  by  Table  3:  During  the  second  and  third  60-minute 
interval,  the  average  gap  is  reduced  by  0.64  percent  and  0.61  per- 
cent, respectively.  In  contrast  to  this  development,  the  decrease  dur- 
ing the  last  three  intervals  is  only  0.25  percent,  0.29  percent,  and 
0.24  percent.  Fig.  6 illustrates  these  findings,  showing  a decreasing 
average  improvement  with  increasing  computing  time.  Therefore, 
we  cautiously  conclude  that  the  objective  function  values  converge. 
Regarding  different  problem  parameters,  no  significant  impact  on  the 
solution  quality  can  be  observed.  Therefore,  the  MS-VNS  seems  to 
perform  with  similar  quality  for  all  problem  settings. 

5.5.  Experiments  with  instances  from  practice 

In  the  third  set  of  experiments,  we  investigated  the  benefits  of  in- 
troducing vehicles  with  multiple  compartments  over  vehicles  with 
a single  compartment  only.  For  these  experiments,  we  generated 
instances  based  on  data  from  practice,  which  we  obtained  from  a 
recycling  company  based  in  the  City  of  Magdeburg,  Germany.  The 
company  is  responsible  for  picking  up  glass  waste  from  containers, 
serving  127  containers  at  44  different  locations.  Up  to  three  different 
containers  are  available  at  each  location,  provided  for  the  collection 
of  colorless,  green  and  brown  glass.  The  trucks,  which  are  used  for 
picking  up  the  glass  waste  at  the  container  locations,  possess  a single 
flexible  wall  that  can  be  introduced  at  predefined  positions  such  that 
the  total  capacity  of  the  truck  can  be  divided  into  at  most  two  com- 
partments. The  position  of  the  wall  can  be  rearranged  for  each  vehicle 
tour. 

The  company  provided  us  with  a list  which  - over  a period  of  sev- 
eral weeks  - documented  for  each  day  which  location  was  visited, 
which  product  type  (glass  color)  was  picked  up,  and  what  the  respec- 
tive supplies  were.  From  this  data  we  were  able  to  estimate  the  so- 
called  fill  rates,  i.e.  the  amount  of  glass  by  which  a container  is  filled 
per  day. 

Over  a planning  horizon  of  25  days  we  then  deterministically  sim- 
ulated the  development  of  the  supplies  for  each  container  at  each  lo- 
cation, assuming  that  a container  is  always  emptied  on  the  last  day 
before  its  capacity  would  be  exceeded.  The  data  (i.e.  the  supplies  for 
each  container  at  each  location)  that  has  been  obtained  by  the  sim- 
ulation for  each  day  has  then  been  considered  as  the  basic  data  for 
the  definition  of  a specific  routing  problem  instance.  These  25  days 
have  been  selected  for  our  numerical  experiments,  giving  rise  to  25 
different  problem  instances.  The  distances  between  the  locations  and 
between  the  depot  and  the  locations  were  determined  by  means  of 
a commercial  navigation  system  and  taken  as  the  costs  of  the  corre- 
sponding edges.  It  was  ensured  that  the  distance  matrix  satisfies  the 
triangle  inequality. 

The  procedure  resulted  in  instances  where  the  number  of  con- 
tainer locations  varies  between  26  and  34  and  the  number  of  sup- 


plies between  34  and  54.  The  number  of  product  types  is  three  for 
all  instances.  The  MS-VNS  approach  was  applied  to  the  25  instances 
allowing  either  one  or  up  to  two  compartments  for  each  vehicle. 
For  both  scenarios,  all  algorithm  settings  were  identical.  As  termi- 
nation criteria,  500n|P|  iterations  without  improvement  and  n|P|/6 
loops  without  improvement  were  used.  The  average  computing  times 
amounted  to  261  seconds  for  the  CVRP-scenario  and  452  seconds  for 
the  MCVRP-scenario. 

Table  4 presents  the  characteristics  of  each  problem  instance  and 
the  corresponding  results  from  the  experiment,  namely  the  number 
of  locations  (n)  and  the  total  number  of  supplies  (s),  the  total  costs 
of  the  solutions  in  case  only  one  compartment  can  be  used  on  each 
vehicle  (tc.cvrp),  the  total  costs  of  the  solutions  if  up  to  two  (flexible- 
sized) compartments  can  be  introduced  (tc.mcvrp)  and  the  improve- 
ment (impr;  in  percent)  which  can  be  obtained  by  the  latter.  Finally, 
the  corresponding  number  of  tours  for  vehicles  with  one  compart- 
ment (tours.cvrp)  and  with  up  to  two  compartments  (tours.mcvrp) 
are  shown. 

The  results  of  this  set  of  experiments  demonstrate  the  economic 
benefits  of  using  vehicles  with  multiple,  flexible-size  compartments 
over  vehicles  with  single  compartments.  For  these  problem  instances, 
the  costs  of  the  tours  could  be  reduced  by  34.8  percent  on  average. 
Furthermore,  the  number  of  tours  (or  vehicles)  necessary  on  aver- 
age to  pick  up  all  supplies  was  reduced  by  one,  from  3.2  tours  to 
2.0  tours. 

6.  Conclusions 

In  this  paper,  a variant  of  the  multi-compartment  vehicle  rout- 
ing problem  was  introduced  and  studied.  This  problem  is  different 
from  the  ones  previously  discussed  in  the  literature.  In  particular, 
it  includes  compartments  of  flexible  sizes,  allowing  for  the  number 
of  compartments  being  smaller  than  the  number  of  product  types, 
which  have  to  be  transported  separately.  Because  of  this  aspect, 
an  additional  question  arises  concerning  the  assignment  of  product 
types  to  vehicles. 

The  main  contributions  of  this  paper  include  (1)  a presentation 
of  this  new  variant  of  the  MCVRP  which  is  inspired  by  a real-world 
application;  (2)  the  development  of  an  exact  and  a heuristic  solution 
procedure  for  this  problem;  and  (3)  extensive  numerical  experiments 
investigating  the  problem  itself  as  well  as  the  performance  of  both 
suggested  algorithms. 

Concerning  the  exact  approach,  the  results  of  the  experiments 
have  shown  that  only  problem  instances  of  limited  size  can  be  solved 
to  optimality  in  reasonable  computing  time.  Detailed  analysis  of  the 
results  from  these  experiments  provided  insights  into  the  impact  of 
different  problem  parameters:  The  problem  becomes  - not  unexpect- 
edly - harder  to  solve  with  an  increasing  number  of  customer  loca- 
tions, product  types  and  supplies.  However,  it  was  also  shown  that 
the  maximum  number  of  compartments,  which  determines  the  num- 
ber of  product  types  per  vehicles,  has  a significant  impact  on  the 
complexity  of  the  problem.  The  smaller  the  number  of  available  com- 
partments is,  the  harder  the  problem  becomes  to  solve.  This  result 
particularly  emphasizes  the  relevance  of  the  additional  assignment 
decision  encountered  in  this  problem,  namely  the  assignment  of 
product  types  to  be  transported  by  each  vehicle. 

Concerning  the  heuristic  approach,  we  were  able  to  show  that 
it  performs  well  - producing  optimal  or  near-optimal  solutions  for 
small  problem  instances.  For  large  instances  with  up  to  450  demands, 
the  heuristic  finds  presumably  good  quality  solutions  in  reasonable 
time. 

In  a third  experiment,  based  on  instances  from  practice,  the  ben- 
efits of  using  vehicles  with  multiple  compartments  of  flexible  sizes 
over  using  vehicles  with  a single  compartment  were  investigated. 
The  results  for  this  specific  application  have  shown  that  - if  multiple, 
flexible-size  compartments  can  be  introduced  - the  total  costs  of  the 
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Table  4 

Results  for  instances  from  practice. 


Instance 

Parameters 

Total  costs 

Number  of  tours 

n 

s 

tc.cvrp 

tc.mcvrp 

impr  (percent) 

tours.cvrp 

tours.mcvrp 

1 

32 

54 

123.83 

83.375 

32.7 

3 

2 

2 

26 

38 

117.59 

78.054 

33.6 

3 

2 

3 

33 

44 

130.18 

91.188 

30.0 

3 

2 

4 

33 

46 

136.24 

95.543 

29.9 

3 

2 

5 

31 

48 

149.20 

79.975 

46.4 

4 

2 

6 

29 

34 

115.05 

77.243 

32.9 

3 

2 

7 

34 

54 

163.42 

94.532 

42.2 

4 

2 

8 

28 

40 

132.67 

92.671 

30.1 

3 

2 

9 

30 

44 

117.10 

78.699 

32.8 

3 

2 

10 

31 

44 

122.70 

82.282 

32.9 

3 

2 

11 

34 

48 

160.17 

90.26 

43.6 

4 

2 

12 

30 

36 

125.60 

85.696 

31.8 

3 

2 

13 

32 

54 

123.83 

83.475 

32.6 

3 

2 

14 

26 

38 

117.59 

78.054 

33.6 

3 

2 

15 

33 

44 

130.18 

91.188 

30.0 

3 

2 

16 

33 

46 

136.24 

95.543 

29.9 

3 

2 

17 

31 

48 

149.20 

79.975 

46.4 

4 

2 

18 

29 

34 

115.05 

77.243 

32.9 

3 

2 

19 

34 

54 

163.42 

94.532 

42.2 

4 

2 

20 

28 

40 

132.67 

92.671 

30.1 

3 

2 

21 

30 

44 

117.10 

78.699 

32.8 

3 

2 

22 

31 

44 

122.70 

82.282 

32.9 

3 

2 

23 

34 

48 

160.17 

90.26 

43.6 

4 

2 

24 

30 

36 

125.60 

85.696 

31.8 

3 

2 

25 

32 

54 

123.83 

83.475 

32.6 

3 

2 

Minimum 

26 

34 

115.1 

77.2 

29.9 

3 

2 

Average 

31.0 

44.6 

132.5 

85.7 

34.8 

3.2 

2.0 

Maximum 

34 

54 

163.4 

95.5 

46.4 

4 

2 

tours  necessary  for  collecting  all  supplies  can  be  reduced  drastically. 
These  results  indicate  the  necessity  of  dealing  with  and  developing 
effective  solutions  approaches  for  multi-compartment  vehicle  rout- 
ing problems. 
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